Non-classical nucleation in supercooled nickel 
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The dynamics of homogeneous nucleation and growth of crystalline nickel from the super-cooled 
melt is examined during rapid quenching using molecular dynamics and a modified embedded atom 
method potential. The character of the critical nuclei of the crystallization transition is examined 
using common neighbor analysis and visualization. At nucleation the saddle point droplet consists of 
randomly stacked planar structures with an in plane triangular order. These results are consistent 
with previous theoretical results that predict that the nucleation process in some metals is non- 
classical due to the presence of long-range forces and a spinodal. 

PACS numbers: 



INTRODUCTION 

Solidification of liquids usually begins by the formation 
of small clusters or nuclei. Understanding the structure 
of nuclei and the dynamics of nucleation and growth is 
of great scientific and technical interest. Nucleation is a 
highly non-linear process and the interplay between the 
non-linearity and a change in symmetry is not under- 
stood in, for example, the crystallization of super-cooled 
liquids. In addition to the interesting scientific questions 
such as the structure of the critical or nucleating droplet 
and the nature of the initial growth, the form of the nu- 
cleation process can have a profound effect on the mor- 
phology of the final state. Classical nucleation theories^ 
predict that the interior of the critical droplet will have 
the same structure as the stable phase. However Klein 
and Leyvraz(KL)j3, Q have argued that the nucleating 
droplet in supercooled liquids with long range potentials 
is non-classical and can have a structure which is signifi- 
cantly different than that of the stable phase due to the 
presence of a spinodal. They suggest that metals with ef- 
fective long range interactions due to elastic forces would 
be a candidate for exhibiting this process which takes 
place arbitrarily close to a spinodal. The non-classical 
droplet is predicted to have a density close to that of the 
supercooled liquid and a triangular structure in 2 dimen- 
sions. In 3 dimensions either randomly stacked triangular 
planes or a bec structure is predicted depending on the 
details of the potential. In the initial growth phase of the 
droplet it is predicted that a core will develop that has 
the density, and symmetry of the stable solid phase, or 
of a metastable solid phase. 0, 0] 

Simple models of elastic systems exhibiting an Austen- 



itc to Martensite transition have indicated a non-classical 
nucleation pattern. However, there has been no test of 
the prediction of non-classical nucleation in either exper- 
iment or simulations with realistic metal potentials. Here 
we present a computational study examining the nucle- 
ation and growth in metastable liquid nickel using state- 
of-the-art semi-empirical potentials. Our results indicate 
that the nucleation is non-classical and is consistent with 
nucleation near a spinodal. It is important to stress that 
by nucleating or critical droplet we refer to the saddle 
point structure. Physically this is the droplet that just 
reaches the critical_size. This appears in the theory as a 
saddle point. [~ 



METHOD 

To perform our simulations we use the modified em- 
bedded atom method (MEAM) potential, referred to as 
potential 4 in Ref. [8|. We extended the radial cutoff, 
r cut , from 4.0 to 4.8 Ato conserve energy in the liquid 
phase. This alteration does not affect the solid prop- 
erties. A thorough discussion of the MEAM formalism 
can be found in the literature 0, 0, 0] and will not be 
repeated here. Recently, this particular potential was 
shown to describe the properties of liquid nickel quite 
well [JlJ- The calculated melting point was within 9 % 
of experiment. The calculated viscosity of 3.84 mPa-s at 
1775 K was in agreement with the lower of the experi- 
mental range of values ~4.0-6 mPa-s. 

We perform a molecular dynamics quenching sequence 
using an isobaric-isothermal (NPT-ensemble). The MD 
calculations use a time step of 1 fs. Constant tempera- 
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ture is maintained by using a Nose-Hoover |12( thermo- 
stat. The Parinello-Rahman method for pressure control 
is used maintaining the pressure to within ±100 bar. The 
sequence is begun by equilibrating a 3D periodic cell of 
(e.g.) 1500 atoms in the liquid state at a temperature 
slightly above the melting point of the system for 25 ps. 
The samples are cooled at 10 K increments and held at 
that temperature for the length of time corresponding 
to the desired cooling rate. For example, a 10 12 K/s 
cooling rate would be obtained by holding the sample at 
each temperature step for 10 ps. The quench rates for 
this study ranged from 10 11 K/s to 10 13 K/s. A few cal- 
culations are performed using a larger cell of 2048 and 
12000 atoms. 

Averages of the potential energy and volume are deter- 
mined at each temperature after 1 ps. The only excep- 
tion to this is for the quench rate of 10 13 K/s, where the 
average temperature and volume are determined for the 
last 0.5 ps. From this data, we determine whether the 
specimen is undergoing a liquid-to-crystal or a liquid-to- 
amorphous transition. The crystallization produces an 
inflection on the cooling curve and the inflection point de- 
termines the solidification temperature, T TO . The liquid- 
to-amorphous transition produces a change of slope in 
the cooling curve. The glass transition temperature, T ff , 
is determined from the intersection of straight lines fit to 
the cooling curve above and below T g . 

We use common neighbor analysis (CNA) to iden tify 
the local crystalline environment of each atom 0, Q, 
IT5L fl6| . This method provides an efficient way to iden- 
tify the local environment of each atom in the system. 
The results are quite insensitive to small displacements 
of the atoms. Since our CNA is not configured to iden- 
tify planar structures we also perform a visual analysis of 
possible planar structures both surrounding the 3D Bra- 
vais lattice identified by the CNA and free standing(i.e. 
without the Bravais core). 

RESULTS AND DISCUSSION 

In Fig. ^ we present the average atomic energies and 
atomic volumes as a function of temperature for various 
quench rates. For quench rates below 5 x 10 12 K/s a 
rapid decrease in energy and volume is seen in the re- 
gion of temperature between 700-900 K. The equilibrium 
freezing temperature is 1540K|llj. The decrease occurs 
over a smaller range of temperature as the quench rate is 
lowered. In contrast, for quench rates of 5 x 10 12 K/s and 
higher (not shown) only a subtle change of slope in the 
energy or volume is seen. We attribute the rapid decrease 
in energy and volume to crystallization and the change 
of slope to solidification into an amorphous state. Below 
the transition, the curves for both energy and volume 
are parallel, but offset from each other. Thus the specific 
heat and thermal expansion of the crystallized material 



are independent of quench rate (which determines local 
structure) but the actual energy and volume differ pri- 
marily due to quenched defects (vacancies) in the solid. 
Above the transition region note that the curves all lie 
on top of each other. The temperature derivative of the 
liquid energy (specific heat) and volume (thermal expan- 
sion) are independent of quench rate. Hence the local 
structure in the liquid is not affected by the quench rate. 

To determine whether there are any important finite 
size effects we compare the results in systems containing 
12000, 2048 and 1500 atoms. The resulting energy curves 
are shown in Fig. [5] The point where the slope begins to 
increase as the temperature decreases indicates the onset 
of crystallization. Nucleation occurs at a slightly higher 
temperature for a given quench rate as the system size 
increases. This implies that the system size dependence 
for nucleation times is as expected |l| and finite size effects 
are not evident. 

Analysis of the resulting microstructure showed that 
for some conditions the small systems are polycrystalline 
and others were single crystals. Similarly the large sys- 
tems exhibit both polycrystalline and single crystalline 
states. The distribution of the polycrystalline domains is 
similar in all sizes examined. In addition the CNA as well 
as visual inspection(to be discussed below) show similar 
patterns of nucleation. We thus conclude that system size 
plays no important role in the resulting microstructure. 

Using CNA and visual inspection we analyze the nu- 
cleation process in supercooled liquid Ni. Upon cool- 
ing for both the large system (12,000 atoms) and the 
small system (1500 atoms), a cluster of atoms having 
the fee structure and the hep structure appear to nucle- 
ate and grow at approximately t a — 1.36 nanoseconds 
into the quench. We use appear since we are looking for 
the saddle point object calculated in the theory. After 
some time into the growth of this cluster the fee particles 
dominate the structure. One implication of this result is 
that the argument of Alexander and McTague^, does 
not apply to this particular systemQ. This argument 
states that in systems where the densities of the liquid 
and solid are similar the solid phase has a bec structure. 
In addition the prediction of icosahedral structures in 
supercooled liquids as suggested by Frank[18| (Lennard- 
Jones), Jonsson et a/.[l|| (Lennard- Jones, LJ), and Ram- 
prasad et al. U3] (Embedded Atom Method, EAM) is not 
observed. We attribute these differences to the spherical 
nature of the LJ and EAM potentials, which do not in- 
clude angular dependence included in the MEAM poten- 
tial. This characteristic appears to prevent the formation 
of icosahedral ordering in pure nickel. 

As discussed above, the KL theory 0, H| predicts 
stacked triangular planes or bec for the droplet struc- 
ture. The KL theory also predicts the development of a 
classical fec/hep core during the growth phase. Therefore 
we look visually for a halo surrounding the fec/hep struc- 
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FIG. 1: Figure a) and b) represent the atomic energy and 
atomic volume as a function of temperature for cooling rates 
between 10 11 and 5 x 10 12 K/s, respectively. 



FIG. 2: Comparison between a large system (12000 atoms), 
a medium system (2048 atoms) and a small system (1500 
atoms) of the atomic energy as a function of temperature. 
The cooling rates shown are 5 x 10 11 and 10 12 K/s. The 
curves are for a single run. 




ture described above. The possibility of a bee halo was 
excluded by the CNA analysis so we search for stacked 
planes with a triangular structure. These were found 
as illustrated in Fig. [3] We then visually examine the 
system at 18 ps prior to the configuration at t a to inves- 
tigate the possibility that there is a droplet consisting of 
only stacked triangular planes with no fec/hep core. This 
configuration is also found and is illustrated in Fig. 




FIG. 3: A 3D representation of the object found at t a . The 
droplet consists of an fec/hep core surrounded by a halo con- 
sisting of planes with an in plane triangular structure. There 
are 80 atoms in the droplet. The red(green) ball represents an 
fcc(hcp) atom while the blue balls are in a planar structure. 
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FIG. 4: A 2D representation of one plane in the saddle point 
droplet. The droplet consists of 4 randomly stacked planes 
with an in plane triangular order. There are 60 atoms in the 
droplet 



These droplets contain ~ 60 atoms. Finally, we investi- 
gate which of these structures; stacked triangular planes 
with, or without, the core correspond to the saddle point 
droplet predicted by the theory. 0, This is done via 
an intervention technique. We take the system at a time 
which we consider it possible that a nucleation droplet ex- 
ists in the sense of a saddle point solution. We then make 
N copies of the system, scramble the velocities keeping 
the temperature fixed, and then restart each copy. In the 
N — » oo limit a saddle point implies that half of the copies 
will continue to grow and half will disappear. |2C| |2lJ If 
the droplet is in the growth phase then N' >> N/2 
droplets will grow. If the copies are made prior to the 
saddle point droplet then N" >> N/2 droplets will dis- 
appear. Using the configuration at t a = 1.36 nanoseconds 
we find that in all of the copies the droplet grows. We 
also find that the same result is obtained for configura- 
tions chosen at 1.5 and 11.5 ps after t Q . Copies made at 
28.5 ps prior to t Q all contain droplets that disappear. 
However, for configurations chosen at 8.5 and 18.5 ps 
prior to t D half of the droplets in the copies disappear 
and half grow. Since we used N = 10 some spread in 
the location of the saddle point droplet was expected. 
In addition, the saddle point flattens at the top near a 
spinodal. Q Even though there was some spread in the 
location of the saddle point droplet it is important to 
note that the droplets at these two times consist only 
of randomly stacked planes with an in plane triangular 
structure and do not contain an fcc/hcp core. 

In conclusion, our study shows that the local struc- 
ture of the initial (saddle point) nucleation droplet con- 
sists of randomly stacked planes with a triangular or- 
dering. During the initial growth phase a dense core 
develops which has both fee and hep atoms. Further 
into the growth phase the core becomes dominated by fee 



atoms. This nucleation sequence is in direct contradic- 
tion of the argument of Alexander and McTague U 
and appears to exhibit the non-classical behavior that 
Klein and Leyvraz0,0] predicted. This result has impor- 
tant implications for the calculation of nucleation rates 
and the understanding of the morphology in either the 
stable or metastable state that has been formed via nucle- 
ation. It also gives additional evidence that the spinodal, 
or limit of metastability of super-cooled liquids may play 
an important role in determining the local structure of a 
critical droplet. Finally, non-classical droplets have been 
found in Lennard-Jones fluids |22l and colloids. [24^ 
however this is the first evidence for non-classical nucle- 
ation influenced by a spinodal in a model of a molecular 
material such as a metal. 
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